Nature of electron Zitterbewegung in crystalline solids 



o 
o 

(N 

CD 



I 



I 

o 
o 



> 

p 

0^ 
O 
On 
O 



X 



Wlodek Zawadzki* and Tomasz M. Rusin -jE 
* Institute of Physics, Polish Academy of Sciences, Al. Lotmkow 32/46, 02-688 Warsaw, Poland; 
\PTK Centertel Sp. z o.o., ul. Skierniewicka lOA, 01-230 Warsaw, Poland 
(Dated: September 2, 2009) 

We demonstrate both classically and quantum mechanically that the Zitterbewegung (ZB, the 
trembling motion) of electrons in crystalline solids is nothing else, but oscillations of velocity assuring 
the energy conservation when the electron moves in a periodic potential. This means that the nature 
of electron ZB in a solid is completely different from that of relativistic electrons in a vacuum, as 
proposed by Schrodinger. Still, we show that the two-band k.p model of electronic band structure, 
formally similar to the Dirac equation for electrons in a vacuum, gives a very good description of 
ZB in solids. Our results indicate unambiguously that the trembling motion of electrons in solids 
should be observable. 

PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 



The phenomenon of electron Zitterbewegung (ZB, the 
trembling motion) was predicted by Schrodinger in 1930 
as a consequence of the Dirac equation for a free rela- 
tivistic electron Schrodinger observed that the elec- 
tron velocity components, given in the Dirac formalism 
by 4 X 4 number operators, do not commute with the 
Dirac Hamiltonian, so the electron velocity in a vacuum 
is not a constant of the motion also in absence of ex- 
ternal fields. The observability of ZB for electrons in a 
vacuum was debated ever since (see e.g. 0,3). Exper- 
imental difficulties to observe the ZB in a vacuum are 
great because the predicted frequency of the trembling 
is very high: hujz — 2moc^ ~ IMeV, and its amplitude 
is very small: Ac = h/rriQC = 3.86 x 10^"^ A. It was later 
suggested that a phenomenon analogous to ZB should 
exist for electrons in semiconductors if they can be de- 
scribed by a two-band model of band structure [3, [E] ■ In 
particular, an analogy between the behavior of free rela- 
tivistic electrons in a vacuum and that of non-relativistic 
electrons in narrow gap semiconductors (NGS) was used 
to predict that the ZB of electrons in NGS should have 
the frequency tihJz — Eg (where Eg is the energy gap), 
and the amplitude \z = h/m*u, where m* is the effec- 
tive mass and u = (£'g/2m*)^/^ ~ 10*cm/s is a max- 
imum electron velocity according to the two-band k.p 
model 0] . This results in much more advantageous char- 
acteristics of ZB as compared to a vacuum; in particular 
Xz ~ 64A for InSb, 37A for InAs, and 13A for GaAs. 
After the papers of Zawadzki 0] and Schliemann et al. 
Q the ZB of electrons in crystalline solids and other sys- 
tems became a subject of intensive theoretical studies, 
see . A classical phenomenon analogous to the ZB was 
recently observed in macroscopic sonic crystals [loj . 

The physical origin of ZB remained mysterious. As to 
electrons in a vacuum, it was recognized that, being of the 
quantum nature, the phenomenon goes beyond Newton's 
First Law. Also, it was remarked that the ZB is due to 
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an interference of states corresponding to positive and 
negative electron energies [ll[ . Since the ZB in solids was 
treated by a two-band Hamiltonian similar to the Dirac 
equation, its interpretation was also similar. This did not 
explain its origin, it only provided a way to describe it. 
However, it was clear that, since the energy bands result 
from electron motion in a crystalline periodic potential, 
in the final count it is this potential that is responsible 
for the ZB. Our paper treats the fundamentals of electron 
propagation in a periodic potential and elucidates the 
nature of Zitterbewegung in solids. The physical origin of 
ZB is of great importance because it resolves the essential 
question of its observability. The second purpose of our 
work is to decide whether the two-band k.p model of 
the band structure, used until now to describe the ZB in 
solids, is adequate. 

It is often stated that an electron moving in a peri- 
odic potential behaves like a free particle characterized 
by an effective mass to*. The above picture suggests 
that, if there are no external forces, the electron moves 
in a crystal with a constant velocity. This, however, 
is clearly untrue because the electron velocity operator 
i^i — Pi/iT^o does not commute with the Hamiltonian 
H = p^/2too -|- V{r), so that Vi is not a constant of the 
motion. In reality, as the electron moves in a periodic 
potential, it accelerates or slows down keeping its total 
energy constant. This situation is analogous to that of a 
roller-coaster: as it goes down losing its potential energy, 
its velocity (i.e. its kinetic energy) increases, and when 
it goes up its velocity decreases. We demonstrate below 
that the electron velocity oscillations due to its motion in 
a periodic potential of a solid are in fact the Zitterbewe- 
gung. Thus the electron Zitterbewegung in solids is not 
an exotic obscure phenomenon — in reality it describes 
the basic electron propagation in periodic potentials. 

The first argument relates to the trembling frequency 
LOz- The latter is easy to determine if we assume, in 
the first approximation, that the electron moves with a 
constant average velocity v and the period of the po- 
tential is a, so LOz — li^vja. Putting typical values 
for GaAs: a = 6.4A, -D = 2.5 x lO^cm/s, one obtains 
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JvjJz = 1.62eV, i.e. the interband frequency correspond- 
ing to the energy gap Eg ~ 1.5eV. The interband fre- 
quency is in fact typical for the ZB in sohds. Next 
we describe the velocity oscillations classically, assuming 
for simplicity a one-dimensional periodic potential of the 
form V{z) = Vq sin(27rz/a). The first integral of the mo- 
tion expressing the total energy is: E ~ mQv'^/2 + V{z). 
Thus the velocity is 
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One can now separate the variables and integrate each 
side in the standard way. However, trying to obtain an 
analytical result we assume V{z) ~ E/2, expand the 
square root retaining first two terms and solve the re- 
maining equation by iteration taking in the first step a 
constant velocity Vzq = (2£'/mo)^/^. This gives z = v^ot 
and 
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Thus, as a result of the motion in a periodic potential the 
electron velocity oscillates with the expected frequency 
ujz = ^TTVzo/a around the average value Vzo- Integrat- 
ing with respect to time to get an amplitude of ZB we 
obtain Az = Voa/{4:nE). Taking again Vq ~ E/2, and 
estimating the lattice constant to be a ~ hpcv/{moEg) 
(see Luttinger and Kohn [IB]), we have finally Az ~ 
fiPcv/ (SiTTmoEg), where Pcv is the interband matrix ele- 
ment of momentum. This should be compared with an 
estimation obtained previously from the two-band k.p 
model 0: Az ~ = h/m*u = h{2/m*EgY/^ ~ 
2hpcv/'moEg. Thus the classical and quantum results de- 
pend the same way on the fundamental parameters, al- 
though the classical approach makes no use of the energy 
band structure. We conclude that the Zitterbewegung in 
solids is simply due to the electron velocity oscillations 
assuring the energy conservation during motion in a pe- 
riodic potential. 

Now we describe the electron velocity oscillations us- 
ing a quantum approach. We begin with the periodic 
Hamiltonian H — p^/2mo + V{z). The velocity operator 
is Vz — Pz/mo- Using the above Hamiltonian one obtains 
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dvz 
dt 
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which is a quantum analogue of the Newton law of motion 
in an operator form. In order to integrate Eq. ([3]) we 
assume a particularly simple periodic saw-like potential. 
It is described by V{z) = —gz for < z < a/2, a < 
z < 3a/2, etc., and V{z) = —Vq + gz for a/2 < z < 
a, (3a/2) < z < 2a, etc., where 17 is a constant force, 
see Fig la. In each half-period z is counted from zero. 
The derivatives are —dV/dz — g ior < z < a/2, a < 
z < 3a/2, etc., and —dV/dz — —g for a/2 < z < a, 
(3a/2) < z < 2a, etc., as illustrated in Fig. lb. Thus 
the electron moves initially with a constant acceleration 
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FIG. 1: Plot of (a) potential, (b) acceleration, (c) velocity, 
and d) position for the electron motion in a saw-like potential 
(schematically) . 



g/ruQ from z = to z = a/2, reaches the maximum 
velocity Vm = {ag/m.oy^'^ at a time tm = {"moa/ gY^^ , 
and then slows down reaching ?) = at a time 2tm- Then 
the cycle is periodically repeated. We calculate v{t) = 
{g/mo)t and z{t) = {g/2ino)t^ in the first, third, fifth 
time (or distance) intervals, and v{t) = Vm — {g/'rno)t 
and z{t) — Vmt — ((?/2mo)t^ in the second, fourth, sixth 
time (or distance) intervals, as illustrated in Figs. Ic and 
Id. We assumed for simplicity v{0) = z(0) = 0. 

It is seen from Fig. Ic that the velocity oscillates in 
time around the average value Vav = Vm/^- This is re- 
flected in the oscillations of position Az(t) around the 
straight line z — Vavt, as seen in Fig. Id. The latter 
are easily shown to be Az = ±a/8, which compares well 
with the classical results. Using the above estimation for 
a we identify again the motion due to periodicity of the 
potential with the Zitterbewegung calculated previously 
with the use of two-band k.p model. 

Strictly speaking, the physical sense of the above op- 
erator reasoning is reached when one calculates average 
values. Our procedure follows the original approach of 
Schrodinger ll|, who integrated operator equations for 
v{t) and z{t). Similar approach is commonly used for 
z(t) and p{t) operators in the harmonic oscillator prob- 
lem [l3|. Equation of motion Q does not contain the 
total electron energy but it is clear from Fig. la that, 
in a quantum treatment, the electron can move either in 
allowed energy bands below the potential tops at i? < 
or above the tops at E > 0. 

Now we describe the ZB using a rigorous quantum ap- 
proach. We employ the Kronig-Penney delta-like poten- 
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tial since it allows us to calculate explicitly the eigenen- 
ergies and eigenfunctions [H, [U . In the extended zone 
scheme the Bloch function is tpk{z) = e^'^^Ak{z), where 

Ak{z) = e-'^^'Ck {e'"' sin[/3fcz] + sin[/3fe(a - z)]} , (4) 

in which k is the wave vector, Cfc is a normalizing con- 
stant and (3k — \/2rnoE /h is a solution of the equation 

sin(/3fca 



cos(/3fca) = cos(fca), 



(5) 



with Z > being an effective strength of the potential. In 
the extended zone scheme, the energies E(k) are discon- 
tinuous functions for k — m:/a, where n = . . . — 1, 0, 1 . . .. 
In this convention, iinn/a < k < {n+l)TT/a, the energies 
E{k) belong to the n-th band and the Bloch states are 
characterized by one quantum number k. Because Ak{z) 
is a periodic function, one may expand it in the Fourier 
series Ak{z) = ^„ A„ exp(ifc„2;), where fc„ = 2nn/a. 

In the Heisenberg picture the time-dependent velocity 
averaged over a wave packet f{z) is 



m) = — I Idkdk' {f\k){k\^\k'){k'\f)e^(^>'-^^'^'/\ 
mo J J idz 

(6) 

where \k) is the Bloch state. The matrix elements of 
momentum are {k\p\k') — h5k',k+k„K{k,k'), where 

K{k,k')^ f\,{zy-^M^dz, (7) 
Jo loz 

The wave packet f{z) is taken in a Gaussian form of the 
width d and centered at /cq, and its matrix elements are 
=E„^ni^(fc:fc»), where 



F(k,kr, 



(8) 



Inserting the above matrix elements to Eq. ^ we obtain 

— V ffdkdk'A*A„'F*(k,k„)F(k',k„')x 
mo ,,JJ 

K{k,k')e^(^^-^y'>'/Hk'M+k,.i9) 



m) 



Figure 2 shows results for the electron ZB, as computed 
for a superlattice. The electron velocity and position 
are indicated. As follows from the inset in Fig. 3, a 
relative narrowness of the wave packet in k space cuts 
down contributions from k values away from ko , and one 
deals effectively with the vicinity of one energy gap. It 
is seen that for a superlattice with the period a — 200A 
the ZB displacement is about ±50A, i.e. a fraction of 
the period, in agrement with the rough estimations given 
above. The period of oscillations is of the order of several 
picoseconds. 

Finally, we want to demonstrate that the two-band 
k.p model, used until present to calculate the Zitter- 
bewegung Q, is adequate for a description of this phe- 
nomenon. We calculate the packet velocity near the point 
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FIG. 2: Calculated electron ZB velocities and position in a 
superlattice versus time. The packet width is d = 400A, 
Kronig-Penney parameter is Z — I.Stt, superlattice period 
is a = 200A. (a) Packet centered at ko = vr/a; (b) and (c) 
packet centered at ko = 0.757r/a. The dashed lines indicate 
motions with average velocities. 
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FIG. 3: ZB of electron velocity in a periodic lattice versus 
time. Solid line: the complete Kronig-Penney model, dashed 
line: the two-band k.p model. Inset: Calculated bands for 
the Kronig-Penney (solid line) and the two-level k.p model 
(dashed line) in the vicinity of fc = -k /a. The wave packet f{k) 
centered at fco = 0.757r/a is also indicated (not normalized). 
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ko = ir/a for a one-dimensional Kronig-Penney periodic 
Hamiltonian using the Luttinger-Kohn (LK) representa- 
tion The LK functions Xnkiz) = e"^^- Unko{z) , where 
Unka{z) = 'u„fe„(z-f a), also form a complete set. We have 



m) = uwrn = 



kk'nn 



{f\nk){nk\v{t)\n'k'){ n'k'\f), 

(10) 

where the velocity in the Heisenberg picture is v{t) = 
(?i/mo)e*^'/''(9/z9z)e~'^*/''. Restricting the above sum- 
mation to the conduction and valence bands: n, n' = 1, 2, 
we obtain in the matrix form 



v{t)ii v{t)i2 
vit)2l v{t)22 



(11) 



where fi = (nfc|/), and v{t)nn' sse the matrix elements 
of the time-dependent velocity operator between the LK 
functions. Equation (jlip looks like the k.p approach to 
ZB used previously. 

The k.p Hamiltonian is obtained from the initial pe- 
riodic Hamiltonian in the standard way [l6| . In the 
two-band model the result is 



, fi^k"^ / 21710 + El fikp2i/mo 
''P " I hkpi2/mo h^P/2mQ + E2 



(12) 



where pi2 — P21 are the interband elements of momen- 
tum, and El and E2 are the band-edge energies. The 
velocity matrix at t = is Vkp = dHkp/hdk. The calcu- 
lation of velocity in the Heisenberg picture is described 
inRef. 

In Fig. 3 we compare the ZB oscillations of veloc- 
ity calculated using: (a) real E{k) dispersions resulting 
from the Kronig-Penney model and the corresponding 
Bloch functions of Eq. (H]); (b) two-band E{k) disper- 
sions obtained from Eq. ((T^) and the corresponding LK 
functions. It is seen that, although we take the packet 
not centered at fc = ir/a, the two-band k.p model gives 
an excellent description of ZB. For fco — n/a the two de- 
scriptions are almost identical. It is seen from Fig. 3 that 
for a wave packet of the width d = 2a the ZB is already 
well described by the k.p model. On the other hand, for 
a much wider packet in k space the two descriptions dif- 
fer more, especially when the fc-width encompasses more 
than one energy gap. 



A few remarks are in order. The transient character 
of ZB, illustrated in Figs. 2 and 3, is a result of describ- 
ing the electron dynamics with the use of wave packets, 
see [1^. In particular, it is seen that wider packets (in 
real space) result in longer transient times. In the lim- 
iting cases of plane waves the electron oscillates indefi- 
nitely, similarly to the classical description [7]. Second, 
one should bear in mind that the standard conductiv- 
ity theories use average electron velocities v = hk/m*, 
as indicated in the above figures. Our work shows that 
at very short times the electron dynamics is completely 
different from its average behavior. Third, both periods 
and amplitudes of ZB, as shown in Fig. 2, are compara- 
ble to those appearing in the Bloch oscillation measure- 
ments [TsJ], so the ZB should be also observable experi- 
mentally. Clearly, it is difficult to follow the behavior of 
a single electron and, in order to observe the trembling 
motion, one should produce many electrons moving in 
phase. This can be most readily done using laser pulses, 
see 0,13. 



In summary, we considered fundamentals of electron 
motion in periodic structures and showed that the exten- 
sively studied phenomenon of electron Zitterbewegung in 
crystalline solids is caused by oscillations of velocity as- 
suring the total energy conservation as an electron moves 
in a periodic potential. This means that, although the 
ZB in solids was often studied in literature using the two- 
band k.p model of band structure analogous to the Dirac 
equation for relativistic electrons in a vacuum, the origins 
of ZB in s solid and in a vacuum are completely differ- 
ent. We also performed a rigorous quantum calculation 
of ZB for an electron in the Kronig-Penney potential and 
showed that the two-band k.p model is adequate for its 
description. 



We dedicate this work to the memory of Professor 
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of Crystalline Solids" was very helpful in our endeavor. 
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